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Abstract 

In this note we review the Singular library INVERSE-SYST.LIB that implements 
Macaulay’s correspondence and other related constructions for local rings. 


1 Introduction 

Let k be an arbitrary field. Let R = k[xi,... x n ] be the ring of the formal series with 
maximal ideal m = (xi, ■ • • , x n ) and let S = k[x 1; ..., x n ] be a polynomial ring, we denote 
by m = (xi,..., x n ) the homogeneous maximal ideal of S. 

In 1916 Macaulay stated a one-to-one correspondence between Artin Gorenstein ideals 
of R and polynomials of S, fl3l . This correspondence can be extended to Artin ideals of R 
and finitely sub-f?-modules of S. Recall that Macaulay’s correspondence is a particular case 
of Matlis duality, see Theorem [23] and Proposition 12.91 

Classically Macaulay’s correspondence has been mainly used for studying homogenous 
ideals, ifTII . lH2ll . Recently Macaulay’s correspondence has been applied to the classification 
of local Artin Gorenstein algebras, see [61, |T|, B3, [HI. Most of the examples appearing in 
these papers have been computed by using Singular, [3l. 

In this note we review the main commands of the Singular library INVERSE-SYST.LIB 
that we used for these computations, 0. The main purpose of this library is to implement 
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Macaulay’s correspondence if the action of R in E r{ (k) — S is defined by differentiation or 
contraction, Theorem l2.21 We also implement some useful operations of S as A-modulc. See 
Section®for a listing of all commands of INVERSE-SYST.LIB. In Section® we give a new 
proof of the classification of Artin Gorenstein local rings with Hilbert function {1, 3, 3,1} 
by using INVERSE-SYST.LIB and the Weierstrass equation of an elliptic curve instead of 
Legendre equation as it was done in |[6l. 


2 Macaulay’s correspondence 

Let A = R/1 be an Artin ring with maximal ideal n = m//. The Hilbert function of A 
is the numerical function HF^ : N —* N defined by HFq(/) = dinik(n*/n l+1 ), i > 0. 
The socle degree of A is the last integer s such that HFa(s) ^ 0. The socle of A is the 
k-vector subspace of A soc(A) = (0 :,4 n), and the Cohen-Macaulay type of A is r(A) = 
dim k (soc(A)). 

Definition 2.1. An Artin ring A with socle degree s is level if soc( A) = tv\ In particular, A 
is Gorenstein iffr(A) = 1. 

The function isAG (I) returns —2 if the quotient A = R/I is not Artin, returns —1 
if A is Artin but not Gorenstein and returns the socle degree of A if the ring A is an Artin 
Gorenstein ring. A similar function is implemented for checking if A is level, see Section® 
The functions socle and cmType compute the socle and the Cohen-Macaulay type of A. 


II we define the ring r 
>ring r=0,(x(1..3)),ds; 

// loading the libary 

>LIB "inverse-syst-v.4.lib"; 

// the quotient r/i is not Artin: 

>ideal i=x(1)"2+x(2)"3, x(2)'4; 

>isAG(i); 

-2 

// the quotient r/i is Artin but Gorenstein: 

>ideal i=x(1)"2+x (2)"3, x(2)"4+x(1)~2, x(3)"2+x(1)*x (2), 
x(1) *x(2)"2*x (3) ; 

>isAG(i); 

-1 

// the quotient r/i is Artin Gorenstein and we get 
// the socle degree 

> ideal i=x(1)"2+x(2) "3, x(2)"4+x (1)"2, x (3) "2+x(1)*x(2); 

> isAG(i); 

4 

// we define an Artin no Gorenstein ideal: 

> ideal i=x(1)"2+x (2) "3, x(2)"4+x (1)"2, x(3)"2+x (1)*x (2), 

x(l)*x(2)"2*x(3); 

> isAG(i); 

-1 

// we compute the socle ideal of r/i 
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> socle (i); 

_[l]=x(l)~2 

_[2]=x(1)*x(2)+x(3)"2 
_[3]=x(2) A 3 
_[4]=x(2) A 2*x(3) 

_ [ 5]=x(1)*x(3)“2 

_[6]=x(2)*x (3) "2 

_ [ 7]=x (3)"3 

// we compute the Cohen-Macaulay type of r/i 

> cmType (i); 

3 

The polynomial ring S can be considered as /I'-modulc with two linear structures: by 
derivation and by contraction. If char (it) = 0, the /^-module structure of S by derivation is 
defined by 



where for all a, (3 e N n , a\ = n"=i “*• 

If char( k) > 0, the f?-module structure of S by contraction is defined by: 



In Singular we can use the above products as follows: 

> ideal F=x(1)~2*x (3)~4+ x(2)"3*x(1)*x (3)+x (2)"5; 

> diff(x(l)"2,F); 

_[1,1]=2*x (3)"4 

> contract(x(1)"2,F); 

_[1,l]=x(3)"4 

It is easy to prove that for any field k there is a /i'-modulc homomorphism 


a : ( S , der ) 


—> (S, cont ) 
h -> al x a 


If char( k) = 0 then a is an isomorphism of J?-modules. The f?-module S is the injective 
hull E r { k) of the f?-module k: 

Theorem 2.2. (fJU\D If k is of characteristic zero then 


E r ( k) = ( S,der ) = ( S,co7it ). 


If k is of positive characteristic then 


E r { k) = ( S, cont ). 
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Given a commutative ring R we denote by Rjmod, resp. Rjmod.Noeth, Rjmod.Artin, 
the category of /(-modules, resp. category of Noetherian /’-modules, Artinian /’-modules. 
The Matlis dual of an /(-module M is M v = Hom R (M, Er(\s.)). We write (—) v = 
Horrift(—, Efiik) ), which is an additive contravariant exact functor from the category Rjmod 
into itself. 


Theorem 2.3 (Matlis duality). The functor (—) v is contravariant, additive and exact, 
and defines anti-equivalence between Rjmod.Noeth and Rjmod.Artin (resp. between 
Rjmod.Artin and Rjmod.Noeth). The composition (—) v o (—) v is the identity functor 
of Rjmod.Noeth (resp. Rjmod.Artin). Furthermore, if M is a R-module of finite length 
then C r (M v ) = £ r (M). 

From the previous result we can recover the classical result of Macaulay, |[T3l . for the 
power series ring, see f9l, lfT2ll . If / C R is an ideal, then ( R/I) v is the sub-R-module of S 

= {g e S | log = 0}, 

this is the Macaulay’s inverse system of I. Given a sub-/(-module M of S then dual M v is 
an ideal of R 

M ± = {f eR I fog = 0 for all g G M}. 


Proposition 2.4 (Macaulay’s duality). Let R = k[xi,... x n j be the n-dimensional power 
series ring over a field k. There is a order-reversing bijection _L between the set of finitely 
generated sub-R-submodules of S = k[?/],... y.„] and the set ofm-primary ideals ofR given 
by: if M is a submodule of S then M 1 - = (0 \r M), and I ± = (0 \s I) for an ideal I C R. 

Given a polynomial H G S of degree l we denote by top(//) the degree l homogeneous 
form of H. 


Proposition 2.5 (Proposition 3.7 and Corollary 3.8 IflTTl . If2ll). Let A = R/1 be an Artin ring 
of socle degree s and Cohen-Macaulay t. The following conditions are equivalent: 

(i) A is level, 

(ii) 1 1 is generated by t polynomials II\..... If 6 S such that deg (Hf) = s, for i = 
1,..., t, and the homogeneous forms top^i),..., top(H t ) are k— linear independent. 

In particular, A = R/I is Gorenstein of socle degree s if and only if 1 1 is a cyclic R-module 
generated by a polynomial of degree s. 

Given a collection of polynomials Hi,... ,H t e S we denote by {Hi,..., Hf) the sub- 
/7-module of S generated by Hi,..., H t . Notice that (i/i, ..., Hf) is not an ideal of S, is 
the k-vector space generated by the collection Hi,... ,H t and their derivatives of any order. 
Notice that in Singular ideals are handled by the list of a given system of generators. In 


4 


the library INVERSE-SYST.LIB the sub-/?-modules of S are handled by using the Singular’s 
structure of ideal, i.e. by the list of a given system of generators. 

We denote by S<i (resp. S <t , resp. S, : ), i <E N, the k-vector space of polynomials of S of 
degree less or equal (resp. less, resp. equal to) to i, and we consider the following k-vector 
space 

i T ±^ __ i ± n s<i + s K i 

(■' )i q ■ 


Proposition 2.6. For all i > 0 it holds 


HF^(i) = dim k (/ _L )j. 


In the library INVERSE-SYST.LIB Macaulay’s correspondence has been programmed 
with respect the two /(-module structures of S. i.e. with respect to the differentiation and 
with respect to the contraction. Here we will show how works Macaulay’s correspondence 
with respect the differentiation. Recall that for technical reasons the sub-/(-modules of S are 
handled in INVERSE-SYST.LIB by using the Singular’s structure of ideal. 

The command invSyst computes the inverse system 7 1 C S' of an ideal / of R, the 
command ideal Ann computes the annihilator M 1 - C R of a finitely sub-/?-module M of 
S. 

In the next example we will show that the composition idealAnn o invSyst is the 
identity map on the set of Artin ideals as Proposition 12.41 predicts. 

> ring r=0,(x(1..3)),ds; 

> LIB "inverse-syst-v.4.lib"; 

// we define an ideal i of r, notice that the first generator is 
// a random polynomial with monomials of degree between 2 and 3 
// and random coefficients between -2 and 2. The second generator 
// is a random homogenous form of degree 3 and coefficients 
// between -1 and 1. 

> ideal i=genPol(2,3,2), genPol(3,3,1), x(2)~3+x(1)*x(3)~4, 

x (1) ~2+x(2) ~2*x(3); 

> i; 

i[1]=2*x(1)~2+2*x(2)~2-x(1)*x(3)+2*x(2)*x(3)-x(3) ~2-2*x (1) ~3 
+x(1)~2*x(2)+2*x(1)*x(2)~2-2*x(2)~3-2*x(1)~2*x(3) 

+2*x(1)*x(2)*x(3)+2*x(2)~2*x(3)-2*x(2)*x(3)~2-x(3)~3 
i [2] =-x (1) "2*x (2) -x (2) " 3+x (1) *x (2) *x (3) +x (2) ~2*x (3) +x (1) *x (3) "2 
+x(3)~3 

i [3] =x (2) "3+x(1)*x (3) "4 
i [4] =x (1) "2+x(2)"2*x (3) 

// we compute the inverse system of i: 

> ideal iv=invSyst(i); 

> iv; 

iv[1]=3*x(1)"2+69*x(2)"2-42*x(1)*x(2)*x(3)-3*x(2)~2*x(3) 

-42*x(1)*x(3)~2+15*x(2)*x(3)"2+22*x(3)"3 
iv [2] =2 4*x (1) *x (3) +3*x (1) *x (2) "2 + 6*x (1) *x (3) ~2-2*x (3) "3 
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// Notice that iv is not cyclic so i is not Gorenstein 
// we compute the annihilator of iv. Should be the ideal i. 

> ideal j=idealAnn(iv); 

> j; 

j[1]=7*x(1)~2+x(1)*x(2)*x(3) 

j [2] =14 *x (2) ~2-7*x (1) *x (3) +14 *x (2) *x (3) -7*x (3) ~2 + 

2 8*x(1)*x(2)"2-52*x(1)*x(2)*x(3)+196*x(2)~2*x(3) 

-98*x(2)*x(3) ~2 

j[3]=343*x(1)*x(2)*x(3)-686*x(2)~2*x(3)+343*x(2)*x(3)~2 
j[4]=5*x(1)*x(3)"2-8*x(2)*x(3)"2+5*x(3)~3 
j[5] =x(2)*x (3) "3 
j [6]=x(3) *4 

// we check that i and j are the same ideal: 

> eqldeal(i,j); 

1 


In the next example we will show that the composition invSyst o idealAnn is the 
identity map on the set of finitely sub-/?-modulcs of S as Proposition 12.41 predicts. 


// We start with a random polynomial 

> ideal q=genPol (2,3,2); 

> q; 

q[l]=2*x(l)~2-2*x(l)*x(2)+2*x(2)~2+2*x(l)*x(3)-2*x(2)*x(3) 

-x (3) ~2-x (1) "3-2*x (1) ''2*x (2) +2*x (1) *x (2) "2-2*x (2) "3 
—2*x(1)*x(2)*x(3)+x(2)'2*x (3)-x (2)*x(3 )~2 
II we compute the annihilator qa of q. 

// The ideal qa is a Gorenstein ideal. 

> ideal qa=idealAnn(q); 

> qa; 

qa[1]=4*x(1)~2-17*x(1)*x(2)-5*x (2)~2 + 12*x(2)*x (3)+x(1)"3 

qa[2]=2*x(1)*x(2)+2*x(2)~2-8*x(1)*x(3)-3*x(1)"2*x(2) 

qa[3]=x(l)*x(3)-x(3)' , 2+2*x(l)*x(2)*x(3) 

qa [4] =x (2) ~3-6*x (1) *x (2) *x (3) 

qa[5]=x(2)'2*x(3)+x(2)*x(3)"2 

qa[6]=x(3)~ 3 

// We get that qa is a Gorenstein ideal and the socle degree 
// of r/i is three that coincides with the degree of q: 

> isAG(qa); 

3 

// we compute the inverse system of qa 

> ideal q2=invSyst(qa); 

> q2 ; 

q2 [l]=6*x(l) *x(2)-24*x(l) *x(3)+22*x(3) / '2 + 17*x(l) "3 
+ 34*x(1)"2*x (2)-34*x (1)*x (2)~2+34*x (2)"3 
+34*x(1)*x(2)*x(3)-17*x(2)"2*x(3)+17*x(2)*x(3)"2 
// from Macaulay's correspondence q and q2 should coincide: 

> eqModlH(q,q2); 

1 
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3 A case study: Artin Gorenstein rings with Hilbert 
function {1,3,3,1} 

As a corollary of the main result of ( 6 ]| we got the classification of Artin Gorenstein local 
rings with Hilbert function {1, 3, 3,1} by using the Legendre equation of an elliptic curve. 
In this section we recover this classification by using the Weierstrass equation of an elliptic 
curve and the library Inverse-syst.lib. 


Theorem 3.1 ([[ 6 ]|). The classification of Artin Gorenstein local k-algebras with Hilbert 
function 11 F 4 = { 1 , n, n. 1 } is equivalent to the projective classification of the hypersurfaces 
V ( F ) C P£" 1 where F is a degree three non degenerate form in n variables. 

See O for the classification of Artin Gorenstein local k-algebras with Hilbert function 
HFa = {1, 7i, n, 1} for n — 1,2. Assume n = 3. Any plane elliptic cubic curve C C P£ 
is defined, in a suitable system of coordinates, by a Weierstrass’ equation, lfl4l proof of 
Proposition 1.4, 

W a y) = x\x 3 = xf + ax \x\ + bx\ 

with a, h E k such that 4a 3 + 27 b 2 0. The j invariant of C is 


j(a,b) 


1728 


4a 3 

4a 3 + 27 b 2 


It is well known that two plane elliptic cubic curves C % = V(W aubi ) C P£, i — 1, 2, are 
projectively isomorphic if and only if j{a 1 , b\) = j(a 2 , b 2 ). 

For elliptic curves the inverse moduli problem can be done as follows. We denote by 
W ( j ) the following elliptic curves with j as moduli: IF(0) = x\x 3 + x- 2 x\ — xf , W (1728) = 
xfx 3 — xgxl — x 3 , and for j f 0,1728 

W(j) = (j — 1728) (xlx 3 + x\x 2 x 3 — xf) + 36xiXg + xf. 


We will show by using the library Inverse-syst.lib that: 


Proposition 3.2. Let A be an Artin Gorenstein local k-algebra with Hilbert function I IF.i = 
{1, 3, 3,1}. Then A is isomorphic to one and only one of the following quotients of R = 

k[[xi,x 2 ,x 3 ]] ; 


Model for A = R/1 

Inverse system F 

Geometry ofC = V (F) C P£ 


X1X2X3 

Three independent lines 

{x{,xix 3 , x 3 x'i, x%, x'i + X 1 X 2 ) 

x 2 (xix 2 - x^) 

Conic and a tangent line 

(x'l, x' 2 , x 3 + 6 x 1 X 2 ) 

* 3 (XlX 2 - X^) 

Conic and a non-tangent line 

(x‘j,XiX 2 ,xf + x '2 — 3X1X3) 

X2X3 - Xi(xi + x 3 ) 

Irreducible nodal cubic 

(Xg,XiX 2 ,XiX 3 ,X 2 ,Xi + 3X^X3) 

X|X3 — xf 

Irreducible cuspidal cubic 

(Xg,Xi + 3X2X3, XiX3,X2 ^ X2X3 + Xg,XiX2) 

W( 0) 

Elliptic curve j = 0 

(X2 + X1X3, X1X2, Xg — 3 Xg) 

W(1728) 

Elliptic curve j = 1728 

I(j) = (x 2 (x 2 - 2 xi ),Hj,Gj) 

W{j) 

Elliptic cun’e with j f 0,1728 
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with: 

Hj = 6jxix 2 — 144 (j — 1728)a;ia;3 + 72 (j — 1728)x 2 a;3 — (j — 1728) 2 ;r§, and 
Gj = jx\ - 12 (j - 1728)xix 3 + 6 (j - 1728)x 2 a: 3 + 144 (j - 1728):r 2 ; 

1(h) = 1(h) if and only if j i = j 2 . 


The first 7 models can be obtained from the corresponding inverse system F by using the 
command idealAnn. Assume that j 7 ^ 0,1728. Let J(j ) be the ideal (W(j)) L \ a simple 
computation shows that = {1, 3, 3,1}, Proposition [276} 


// we define a ring of characteristic zero, three variables and ground field 
// a field of functions with indeterminate c(l) 

> def r=workringc(0,1,3); 

> setring r; 


> r; 





// 

characteristic 

0 


// 

1 parameter 


c (1) 


// 

minpoly 


0 


// 

number of vars 

3 


// 

block 

1 

ordering 

ds 

// 



names 

x ( 

// 

block 

2 

ordering 

C 


// the ideal p defines a elliptic curve with j=c(l) 

> ideal p=weierstrassp (); 

> p; 

p [1] = (-c (1) +1728) *x (1) ''3+ (c (1) -1728) *x (1) *x (2) *x (3) + 

(c(1)-1728)*x(2)~2*x(3)+36*x(l)*x(3)~2+x(3)"3 
// we define the ideal q. We will prove that the inverse system of q 
// is p. 

> ideal q=idealwp(); 

> q; 

q [ 1 ] = (6*c (1) ) *x (1) *x (2) + (-144*c (1) +2 48832) *x (1) *x (3) + 

(72*c(l)-124416)*x(2)*x(3)+(—c(1)~2+3456*c(l)-2985984)*x(3)~2 
q[2] = (c(l))*x(l) - '2+(-12*c(l)+20736)*x(l)*x(3) + (6*c(l)-10368)*x(2)*x(3) + 
(144*c(l)-248832)*x(3 )~2 
q[3]=-2*x(1)*x(2)+x(2)"2 

// we check that q is contained in p~\perp 

> dif f (q, p) ; 

_[1,1]=0 

_[1,2]=0 
_[1,3]=0 

// If we perform de division of the $4$-th power of the maximal ideal 
// by q we get three matrices Q, R, U such that (see Singular's manual) 

// generators(maxideal(4))*U=generators(q)*Q + R 

// U is the 15x15 identity matrix, R is the 15x1 zero matrix 
// and Q is a 6x15 matrix with coefficients in the ground field (see below 
// for more details). The command is: 

> division(maxideal(4),q); 

From the last computation we get that the denominators of the coefficients of the matrix Q 
are constant polynomials or polynomials with roots in {0,1728}. Hence for all j = c(l) f 


8 



0,1728 we get that m 4 C q, so q is an Artin ideal. Notice that for all j = c(l) G k, 
q = I(j ) and p = (W(j)). Since I(j ) is generated by three homogeneous elements, I (j ) 
is a homogeneous complete intersection ideal. In particular I(j ) is a homogeneous Artin 
Gorenstein ideal, so HFr//^) is symmetric. Notice that the generators of /(j) are three 
homogeneous forms of degree two, so the Hilbert function of A = R/I(j ) is {1, 3, 3,1}. We 
know that I(j) is contained J(j ) = (W (j)) x . Since HF R/j(j) = HF R/j(j) = {1, 3, 3,1}, we 
get that I(j) = J(j) = (W{j)) L , i.e. /(j) = {W(j)) ± . 

4 Commands 

Next, we list the most important commands of Inverse-SYST.lib. 

IDEAL THEORY 
genPol(i,j,a). 

Returns a generic polinomial sum of forms of degrees 
between i and j, with integer coefficients in [-a,a]. 
eqldeal(J,I); I, J ideals. 

Returns 1 if I=J, 0 otherwise, 
socle(J); I ideal 

Returns -1 if J is not Artin, 
returns the ideal of J if J is Artin. 
cmType(J); J ideal. 

Returns -1 if J is not Artin, 

returns the Cohen-Macaulay type of J otherwise. 
isAG(I); I ideal. 

Returns -2 if J is not Artin, 

returns -1 if J is Artin but not Gorenstein, and 
returns the socle degree if J is Artin Gorenstein. 
isLevel(I); I ideal. 

Returns -2 if J is not Artin, 

returns -1 if is Artin but not Level, and 

returns the socle degree if J is Artin Level 

MACAULAY INVERSE SYSTEM CORRESPONDENCE WITH COEFFICIENTS 

invSystG(ideal J) 

returns the inverse system of J; J is Artin Gorenstein 
invSyst(J) 

returns the inverse system of J; J is Artin 
idealAnnG(poly f) 

returns the Artin Gorenstein ideal with inverse system f 
idealAnn(I) 

returns the Artin ideal with inverse system I 

STRUCTURE OF INJECTIVE HULL WITH COEFFICIENTS 

memberIH(g,I); I=fl,... list of polynomials, g polynomial. 

returns 1 if g belongs to the R-submodule of S generated by fl,...fs in S, 
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returns 0, otherwise 

subModlH(I,J); I=fl,... list of polynomials, J=gl,... list of polynomials. 
Returns 1 if I is a sub-R-submodule of J, both sub-R-modules of S; 

0 otherwise 

eqModIH(I,J); I=fl,... list of polynomials, J=gl,... list of polynomials. 
Returns 1 if I=J, both sub-R-modules of S; 

0 otherwise 

minGensIH (I); I=fl,...,fs list of polynomials. 

Returns a minimal system of generators of <fl,...,fs>, sub-R-module of S 
colonlnvSyst(f,g); f, g polynomials. 

Returns an element g of R such that Aof=g if exists, 

0 otherwise. 

MACAULAY INVERSE SYSTEM CORRESPONDENCE WITH NO COEFFICIENTS 
invSystGNC(ideal J) 

returns the inverse system of J; J is Artin Gorenstein 
invSystNC(J) 

returns the inverse system of J; J is Artin 
idealAnnGNC(poly f) 

returns the Artin Gorenstein ideal with inverse system f 
idealAnnNC(I) 

returns the Artin ideal with inverse system I 

STRUCTURE OF INJECTIVE HULL WITH NO COEFFICIENTS 

memberlHNC(g,I); I=fl,... list of polynomials, g polynomial. 

Returns 1 if g belongs to the R-submodule of S generated by fl,..,fs in S, 
returns 0 otherwise 

subModIHNC(I,J); I=fl,... list of polynomials, J=gl,... list of polynomials. 
Returns 1 if I is a sub-R-submodule of J, both sub-R-modules of S; 
returns 0 otherwise 

eqModIHNC(I,J); I=fl,... list of polynomials, J=gl,... list of polynomials. 
Returns 1 if I=J, both sub-R-modules of S; 
returns 0, otherwise 

minGensIHNC(I); I=fl,...,fs list of polynomials. 

Returns a minimal system of generators of <fl,...,fs>, sub-R-module of S 
colonlnvSystNC(f,g); f, g polynomials. 

Returns an element g of R such that Aof=g if exists, 
returns 0 otherwise. 

RINGS WITH PARAMETERS 

workringc(p, t, n) 

returns the def of a ring r with 
t coefficients c (1),...,c (t), 
n vars x(1),...,x (n), char is p, and order ds 

ELLIPTIC CURVES 

weiertrass j (t) 
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returns the ideal generated by Weierstrass equation 
of the elliptic curve with j invariant 
idealwj(t) 

returns the ideal with inverse system weierstrassj(j) 
weiertrassp () 

returns the ideal generated by Weierstrass equation of the 
elliptic curve with moduli j=c(l) 
idealwp(t) 

returns the ideal with inverse system weierstrassp() 
with moduli j=c(l) 
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